Electrically tunable three-dimensional ^-factor anisotropy 
in single InAs self-assembled quantum dots 
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Three-dimensional anisotropy of the Lande ^-factor and its electrical modulation are studied for 
single uncapped InAs self- assembled quantum dots (QDs). The ^-factor is evaluated from measure- 
ment of inelastic cotunneling via Zeeman substates in the QD for various magnetic field directions. 
We find that the value and anisotropy of the ^-factor depends on the type of orbital state which arises 
from the three-dimensional confinement anisotropy of the QD potential. Furthermore, the ^-factor 
and its anisotropy are electrically tuned by a side-gate which modulates the confining potential. 



The Lande ^-factor, the magnetic response of spin, is 
a physical constant reflecting the spin-orbit interaction 
(SOI) and the quantum confinement effect in semicon- 
ductor nanostructures, since it is determined by coupling 
between orbital and spin angular momentum^. In low- 
dimensional systems, the confining potential has large 
asymmetry resulting in an anisotropic ^-factor. There- 
fore, the ^-factor anisotropy can be electrically modu- 
lated by gating the quantum dot (QD) confining po- 
tential. This may be exploited for coherent manipula- 
tion of electron spins through g-tensor modulation reso- 
nance (g-TMR) which was previously studied for a quan- 
tum well 2 . For single electron spins, self- assembled QDs 
(SAQDs) and nanowire quantum dots (NWQDs) made 
out of narrow gap semiconductors such as InAs, InP and 
InSb are relevant for the study of ^-factor anisotropy 3 - - — , 
because they have quite large negative values of ^-factor 
for electrons due to the strong SOI. In particular for InAs 
QDs we previously demonstrated that both the SOI effect 
and the orbital states are influenced by three-dimensional 
(3D) electrostatic potential^. 

Among SAQD systems InAs SAQDs are the most ex- 
tensively studied in crystal growth as well as optical and 
electrical characterization. The InAs SAQDs in our study 
are uncapped or unstrained so that the QD size is rela- 
tively large, laterally 100 nm wide and vertically 30 nm 
high. The QD shape is anisotropic with the confine- 
ment strong in the out-of-plane direction and weak in the 
in-plane direction, leading to ^-factor anisotropy^. Fur- 
thermore, the in-plane confinement is so weak that the 
confinement potential and the confined electron wave- 
function can be modulated by means of electrical gating. 
A local or anisotropic gating has been applied to InAs 
SAQDs£"i£ and NWQDs±± to modulate in- situ QD-lead 
tunneling coupling, and angular anisotropy of SOI en- 
ergy and g-factor. In our previous study-iS of the elec- 
trical tuning of ^-factor we identified tunability of the 
g-tensor only in a two-dimensional (2D) plane and as- 



sumed that the QD could be approximated as a disk-like 
2D harmonic potential as if often don o 12 i 13 . Studies of 
the anisotropy of the SOI have however shown that a 3D 
confinement, arising from the QD shape as well as the 
metal electrodes asymmetrically contacted to the QD, 
may be more realistic 8 . Such a 3D confinement leads to 
arbitral direction of orbital angular momentum (OAM) 
and therefore arbitral ^-factor anisotropy^ which can be 
modulated largely by electrical gating. 

In this work, we investigated full 3D ^-factor 
anisotropy for single InAs SAQDs. The ^-factor is de- 
rived by measurement of inelastic cotunneling or Kondo 
effect through the Zeeman substates in the QD 15 i 16 , since 
^-factor cannot be evaluated adequately via magnetic 
evolution of ground states due to the large electron charg- 
ing energy^. The obtained ^-factor shows a different 
3D anisotropy depending on the charge state. From this 
behavior we distinguish the related orbital type as ei- 
ther that with small OAM type (s-orbital like) or large 
OAM (p-orbital like) confined to a 3D potential rather 
than a 2D potential. We used an electrical sidegate 
(SG) placed nearby the QD to largely modulate the 3D 
QD confining potential and therefore the 3D ^-factor 
anisotropy. We finally characterize the performance of 
g-TMR in the present InAs SAQD using the values of 
^-factor anisotropy obtained here. 

Uncapped InAs SAQDs were grown with Stranski- 
Krastanov mode by molecular beam epitaxy on a [001] 
semi- insulating GaAs substrate. A pair of Ti (3 nm)/Al 
(150 nm) source and drain electrodes separated by a 30 
nm gap was deposited on the surface of a single SAQD 
using electron beam lithography techniques^. The QD 
is positioned at the edge of the nanogap and has a larger 
overlap with the drain lead than the source as shown in 
Fig. 1(a). Differential conductance of electron transport 
through the QD was measured for a small source-drain 
voltage V s d at a bath temperature of 40 mK using con- 
ventional lock-in techniques. To change the number of 
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FIG. 1. (a) Scanning electron micrograph image of the device 
measured. An In As QD is contacted to the source and drain 
electrodes at the edge. Only the sidegate nearest the QD is 
used for the present study. The coordinates are defined as 
the figure, (b) Differential conductance as a function of V 3 d 
and Vb g measured for B z = 1 T and V sg — V. The dashed 
line indicates the center of the Coulomb diamonds which is 
shifted due to leakage of the back-gate to the source lead as 
discussed in the main text, (c) Magnetic evolution of inelastic 
cotunneling peaks and the Kondo anomaly in Region I. Below 
100 mT, superconducting features appear because the Al elec- 
trodes are superconducting. The Kondo zero-bias anomaly is 
only observed when \g z \fiBB z < ksTx- (d) Magnetic evolu- 
tion of inelastic cotunneling steps in Region II. The red traces 
in (c) and (d) show the data at B z = 0.8 T. 



electrons in the QD the backgate (BG) voltage Vb g ap- 
plied to a doped layer beneath the QD was modulated. 
The ^-factor was evaluated from the transport measure- 
ment for various magnetic field angles selected in-situ 
using a vector magnet system. 

Figure 1(b) shows the stability diagram with a per- 
pendicular magnetic field B z = 1 T and a SG voltage 
V sg = V. We observe a series of diamond shaped re- 
gions or Coulomb blockade regions. Note that the QD 
of interest is coupled with another small QD which ex- 
ibits very small tunnel currents and has little effect on 
the transport discussed here^. The center of the dia- 
mond is shifted towards negative V s d as the BG voltage 
Vbg is increased. This shift is caused by current leakage 
from the BG to the source electrode and does not influ- 
ence the transport characteristics of the QD, so we regard 
the dashed line as V s d = V hereafter. In three differ- 
ent Coulomb blockade regions, I, II, and III in Fig. 1(b), 
two peaks for I and two steps for II and III are observed 
on both sides to the dashed line. These peaks or steps 
are due to inelastic cotunneling through the QD, which 



appears when eV s d is consistent with excitation energy 
in the QD^. The onset of inelastic cotunneling is usu- 
ally signified by steps rather than peaks. The peaks are 
only observed when the cotunneling is associated with 
the Kondo process or when the peak separation is not 
very large as compared to the Kondo temperature Tk^- 
This is indeed the case for the present QD. 

We measure the magnetic field evolution of the inelas- 
tic cotunneling peaks to study the inelastic cotunneling 
and the Kondo effect for the three regions. The results 
measured for Region I and II are shown in Fig. 1(c) and 
(d), respectively. The Al electrodes have a supercon- 
ducting critical field B c ~ 100 mT and here we focus 
on the field range above 100 mT where the leads are in 
the normal state. Note that two peaks observed at V s d 
~ ± 0.2 mV below 100 mT are caused by quasi-particle 
tunneling&2£ between the superconducting leads. 

In Fig. 1(c), the two inelastic cotunneling peaks merge 
to a zero-bias peak or Kondo peak as the magnetic field 
is decreased from B z = 1 T down to 500 mT. The zero- 
bias peak is unchanged with B z down to 100 mT. From 
the peak width we evaluate Tk ^ 1 K for Region I. 
The two inelastic cotunneling peaks become separated 
approximately linearly with B z , reflecting Zeeman en- 
ergy \g z \fJ>BB z , where \ib is the Bohr magneton. From 
this behavior we evaluate the ^-factor of \g z \ = 3.5 for 
Region I. The Zeeman energy is 0.1 meV at B z = 500 
mT, comparable to ksTx-, where ks is the Boltzmann 
constant. Note that the Zeeman splitting appears only 
when the electron number in the QD is odd, and therefore 
the SU(2) Kondo effect appears at low magnetic fields. 
For Region II the two inelastic cotunneling steps split by 
the Zeeman effect are observed at \eV s d\ = \Qz\^bB z with 
\g z \ = 5.2 in Fig. 1(d). For Region III, similar inelastic co- 
tunneling steps are observed but with a different ^-factor 
of \g z \ = 3.8 (not shown). The Kondo zero-bias anomaly 
is not clear in Region II and III, probably because Tk is 
lower than the measurement temperature. 

Then we set the magnetic field at B = 1 T directed 
with various angles measured from the -\-z axis in the 
x-z and y-z plane and (j) measured from the +x axis in the 
x-y plane, and measured the inelastic cotunneling peaks 
(or steps) split by A(0, (j)) in V s d- We evaluated the g- 
factor at various sets of and (j) as given by eA//i#£?, and 
derived 3D magnetic angular dependence of the ^-factor 
as shown in Fig. 2 (points). Color indicates the abso- 
lute value of the ^-factor in each magnetic field direction 
represented as a distance from the coordinate origin. In 
this figure the ^-factor takes a maximum in the direction 
approximately (0,</>) = (-45°, 180°) and (135°, 0°) tilted 
from the z-axis by —45° in the x-z plane, and takes a min- 
imum in the direction about (45°, 0°) and (-135°, 180°). 

General model for a 3D anisotropic ^-factor 5 is g(B) = 
\Jg\B\ + g\B\ + g\B\l\B\, where (g u g 2 , gs) are the 
\g\ values in the orthogonal principal axes directions and 
(Bi, £?2, Bs) are the magnetic field components along 
the principal axes. These may be defined as B\ = 
B cos(0 + 0o ) sin(0 + O ), B 2 = B sin(</> + O ) sin(0 + O ) 
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various directions and its 3D surface plot (gray) obtained by FIG. 3. The ^-factor in the cut plane of x-y (a), y-z (b), and 

the theoretical fit to the data using the equation in the main x-z (c) of the 3D polar plot of data in Region I at V sg = 

text. The distance of each plotted point from the origin is V presented in Fig. 2. The red line is the theoretical fit. The 

a measure of the absolute value of the ^-factor which is also dotted lines indicate the axes of the maximum and minimum 

indicated by color. ^-factor directions. 



and Bs = B cos(6 + 6q). Here 0o and Oq define the off- 
set of the ^-factor principal axes from the measurement 
axes. Figure 2 also shows a 3D surface plot of the ^-factor 
in Region I with V sg = V fitted by the above equa- 
tion with the following parameters, (#i, #2, #3, #o 5 0o) = 
(5.2, 2.3, 2.2, -49°, 12°). 

The \g\ values in the cut plane x-y, y-z, and x-z in Fig. 2 
are shown in Figs. 3(a), (b), and (c), respectively. The g- 
factor is anisotropic with magnetic field directions in all 
figures and the anisotropy is significantly different from 
figure to figure. In the x-y cut plane of (a), the ^-factor 
takes a maximum (minimum) of \gmax\ =4.1 (\gmin\ — 
2.4) at (j) = —8° measured from the x-axis. Similar shape 
of anisotropy is observed in the y-z cut plane of (b) with 
the maximal ^-factor at = —11° measured from the 
z-axis. The ^-factor anisotropy is the largest in the x-z 
cut plane of (c) with \g ma x\ = 4.9 and \g m in\ = 2.2 at 
~ ^45° measured from the z-axis, respectively. These 
results differ from our previous study for a QD with ge- 
ometrically symmetric coupling to the leads 7 in which 
the ^-factor was isotropic for in-plane rotation and max- 
imal in the z-direction. The ^-factor in QDs strongly 
depends on the symmetry of the confinement which af- 
fects OAM 14 . The larger the OAM the larger the \g\ 
value. Therefore, the OAM of electrons in Region I is 
the largest when the magnetic field is applied in the x-z 
plane with ~ —45°. The same measurement and theo- 
retical fitting of the ^-factor anisotropy as shown in Fig. 2 
were performed for Region II and III. The results for the 
x-z cut plane are shown in Fig. 4(a), and (b) for Region 
II, and III, respectively. The ^-factor anisotropy in Re- 
gion II is similar to that in Region I, showing a maximum 



at ~ -45° with I 

9max\ — 5.4. For Region III, however, 
the ^-factor is almost isotropic with \g\ =4, which is the 
same for the other planes, indicating a spherical distribu- 
tion of the ^-factor. This difference may be accounted for 
by assuming that the orbital type is p-like (large OAM) 
for Region I and II, whereas it is s-like (small OAM) for 
Region III. 

We assume that the tilting of the ^-factor principal axis 
from the z-axis, observed for Region I and II, is due to the 
3D asymmetry of the QD confining potential. In Fig. 1(a) 
the QD is positioned at the edge of the nanogap and ap- 
proximately one half of the QD is covered by the drain 
electrode metal while a small part of the QD is covered 
by the source electrode metal. In our previous work 7 -, we 
discussed that electrons in the QD beneath the source 
and drain electrode metal are depleted. Therefore, the 
potential shape of the QD is considered as a half pyra- 
mid in the uncovered region by the contact electrodes 18 . 
We used an 8-band k-p theory to calculate the prob- 
ability density of electrons confined by such an exotic 
potential. A p-orbital state confined by a 2D potential 
extends over the x-y plane, having the largest magnetic 
coupling with a magnetic field along the z-axis. In the 
present case, however, a p-like orbital state extends over 
a plane tilted by about +45° from the z-axis 18 and there- 
fore the ^-factor principal axis is tilted by —45° from the 
z-axis 14 , which is consistent with our finding in Fig. 3(c) 
and Fig. 4(a). To further examine the 3D confinement, we 
studied the in-plane magnetic field effect on the Coulomb 
peaks 18 . The in-plane magnetic field evolution of states 
is not linear, indicating that not only Zeeman effect but 
also OAM effect are involved and that OAM may not be 
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FIG. 4. Anisotropic ^-factor for Region 11(a) and 111(b) with 
V sg = V in the x-z plane plotted in the same way as shown 
in Fig.3(c). 
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FIG. 5. Anisotropic ^-factor in the x-z plane measured for 
Region I at V sg = -0.5 V in blue, V in red, and +0.4 V in 
green. As V sg is made positive, the ^-factor get small and the 
^-factor maximal direction gmax approaches the z-axis. The 
optimized direction of magnetic field B op t is selected as Q± 
becomes the largest at V sg — V (explained in the text). 



perpendicular to the in-plane magnetic field in the 3D 
confinement potential. Note that in another sample hav- 
ing a QD symmetrically covered by the source and drain 
electrode metal, the observed ^-factor has a principal axis 
along the z-axis and it is almost isotropic with respect to 
the in-plane magnetic field directions^. 

We now discuss use of the SG to modify the ^-factor 
anisotropy. The electric line of force due to the SG volt- 
age V sg extends through the GaAs substrate with its large 
permittivity, so that application of positive V sg may shift 
the electron wavefunction downwards in the direction of 
the SG. Figure 5 shows the ^-factor anisotropy measured 
for magnetic field in the x-z plane at various V sg of - 
0.5, 0, and +0.4 V in Region I. As V sg is made posi- 



tive, the amplitude of the ^-factor becomes small, and 
the ^-factor maximal direction gmax approaches the z- 
axis. This behavior is qualitatively predicted from our 
simulation^ which shows that the OAM vector of the p- 
like orbital state tends to be tilted toward the z-axis so 
does the g m ax direction as V sg is made positive. When 
the electron wavefunction shifts to the bottom of the QD 
the \g\ value may become small, because Ga atoms are 
diffused from the GaAs substrate into the QD near the 
substrate22i22. 

Finally, we estimate the possible Rabi frequency of 
electron spin by means of g-TMR from the obtained re- 
sults. The anisotropic ^-factor is described by a tensor, 
9ij (hj — x iVi z )- For application of a static magnetic 
field Bq, the spin precession vector CIq is given by 
= ^ZjgijBj. Then suppose to apply microwave (MW) 
on the SG. Because the ^-factor anisotropy or g-tensor 
depends on V sg as shown in Fig. 5, the MW induced a.c. 
voltage gives rise to an a.c. precession vector which has 
a component fl± perpendicular to the initial Ho- The 
Q± is maximal for Bq || B op t in Fig. 5. For example, for 
the Larmor precession frequency of 20 GHz with Q and 
MW induce a.c. voltage of 1 mV, the calculated Rabi 
frequency is 0.88 MHz. This value is similar to that in 
our previous work 10 . 

In summary, we have investigated 3D ^-factor 
anisotropy in single InAs SAQDs by measuring inelas- 
tic cotunneling for three different charge states. We 
observed that the ^-factor principal axis is tilted by 
about 45° from the growth direction, and assigned it to 
the confinement asymmetry created by the QD shape 
as well as the geometry of the source and drain elec- 
trodes. We observed anisotropic and isotropic ^-factor 
depending on the charge state related to the orbital type. 
In addition, we succeeded in electrically tuning the g- 
factor anisotropy via the confinement potential modu- 
lation with V sg . The ^-factor anisotropy can be opti- 
mized toward a manipulation of electron spins by means 
of g-TMR in the QDs by appropriately designing the QD 
shape, confinement potential, and the coupling efficiency 
to the SG. 
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